Chaos in circuit QED: decoherence, localization, and nonclassicality 
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We study the open system dynamics of a circuit QED model operating in the ultrastrong cou- 
pling regime. If the resonator is pumped periodically in time the underlying classical system is 
chaotic, implying build-up of seemingly irregular sub-Planck structures of the Wigner phase space 
distribution. Indeed, the periodically driven Jaynes-Cummings model in the Born-Oppenheimer 
approximation resembles a Duffing oscillator which in the classical limit is a well-known example 
of a chaotic system. Detection of the field quadrature of the output field acts as an efi'ective po- 
sition measurement of the oscillator. We address how such detection affects the quantum chaotic 
evolution in this bi-partite system. We differentiate between single measurement realizations and 
ensembles of repeated measurements. In the former case a measurement/decoherence induced local- 
ization effect is encountered, while in the latter this localization is almost completely absent. This 
is in marked contrast to numerous earlier works discussing the quantum-classical correspondence in 
measured chaotic systems. This lack of a classical correspondence under relatively strong measure- 
ment induced decoherence is attributed the inherent quantum nature of the qubit subsystem and in 
particular the sustainable quantum correlations shared between the qubit and the field. 

PACS numbers: 42.50.Pq, 42.50.Lc, 42.65.Sf 



I. INTRODUCTION 

While chaos is well defined in classical mechanics, the 
degree to which it can manifest itself in quantum me- 
chanics remains a problem of fundamental interest [1]. 
For example, classically chaotic systems exhibit expo- 
nential sensitivity to initial conditions, something that 
is absent from quantum dynamics which is manifestly 
unitary. Moreover, the discrete spectrum of closed quan- 
tum systems implies quasi periodicity which contradicts 
the irregular behaviour seen in classical chaotic systems. 
Theoretical work has indicated that the inclusion of de- 
coherence, on the other hand, can lead a quantum system 
to become exponentially sensitive to initial state fluctu- 
ations [2]. 

On the experimental front, the advent of the laser cool- 
ing of atoms has allowed the realization of a number of 
classic model hamiltonians that display quantum chaotic 
phenomena in a simple setting. Of particular relevance to 
this paper we mention the observation of dynamical local- 
ization of momentum in a kicked rotor [3j , where classical 
(chaotic) diffusion is suppressed by quantum interference, 
and chaos assisted tunneling [H [5] , where tunneling oc- 
curs in phase space between islands of regular motion, 
thereby passing through the classically insurmountable 
dynamical barriers provided by the boundaries between 
integrable and chaotic regions. The kicked top has also 
been realized using atomic spins [6 . This latter exper- 
iment used continuous weak measurement techniques to 
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reconstruct the full quantum state as a function of time 
and hence allow the visualization of quantum dynam- 
ics directly in phase space. Some of these experimental 
setups have also been used to investigate the effects of 
decoherence [3 [8], and, as expected, it was found that in 
each case the system reverted to classical behaviour. 

Driven by the promise of technological applications, 
continuing progress in the manufacture of man made 
and hybrid quantum systems also provides signiflcant op- 
portunities for exploring basic quantum dynamics under 
highly controllable conditions [9 . Within the realm of 
this progress the idea of quantum simulators has become 
reality [lQi-il2j . While quantum simulation of closed sys- 
tems has been the main focus thus far, the extension to 
open systems has been considered [T3^ . In our view these 
new quantum systems naturally lend themselves to the 
study of chaotic quantum dynamics subject to decoher- 
ence. 

In this work we explore quantum chaos in a circuit 
quantum electrodynamics (QED) setup comprised of a su- 
perconducting qubit resonantly coupled to a single mode 
of a high-Q transmission line resonator [14 . The pres- 
ence of the two-level qubit means that this is a man- 
ifestly quantum system in the sense that it does not 
possess a simple classical limit. The parameters of the 
corresponding Hamiltonian can be externally controlled 
to a large extent, and even the strength of the decoher- 
ence can be manipulated by various means. In order 
to achieve chaotic dynamics we consider the case where 
the resonator is externally pumped by a drive oscillator 
with adjustable amplitude and phase. By changing ei- 
ther the qubit-fleld coupling [15 , the pump amplitude or 
the pump modulation frequency, the system can be var- 



ied between being regular or chaotic. In particular, in the 
ultrastrong coupling regime ^ , and at vigourous enough 
pumping, the system becomes chaotic in the sense that 
the semi-classical equations-of-motion corresponding to 
a classical electromagnetic field (but retaining the quan- 
tum nature of the qubit), become chaotic. Furthermore, 
in the Born-Oppenheimer approximation (BOA) where 
the excited state of the qubit is adiabatically eliminated, 
we find that the driven circuit QED system mimics a 
quantum Duffing oscillator [16 . The duffing oscillator is 
an example of a well-studied system whose classical limit 
is known to be chaotic. 

In the chaotic regime we shall see that an initially lo- 
calized phase space distribution of the electromagnetic 
field rapidly broadens and the Wigner distribution builds 
up seemingly irregular sub-Planck structures character- 
istic of chaos. In comparison to how a localized phase 
space distribution evolves under the corresponding semi- 
classical equations-of-motion, the broadening in the fully 
quantum evolution is much more dramatic leading to a 
breakdown of the Ehrenfest theorem. One possible route 
to recover the classical dynamics is to try to localize 
the state by weakly measuring the position quadrature 
of the field. This induces decoherence and non-unitary 
time evolution [17] . One of the main results of our work 
concerns the effect of such decoherence on the present 
chaotic bi-partite system. In general, it has been shown 
that for quantum chaotic systems that possess a clear 
classical limit, decoherence is a possible candidate to ex- 
plain how chaos, in a classical sense, emerges in quantum 
mechanics. In particular, decoherence can prevent rapid 
spreading of an initially localized state and the Ehrenfest 
theorem remains valid for long times [2^. However, when 
the system does not support a classical limit, like here, 
the situation becomes more complex [18]. 

In the present model the effects of quantum ffuctua- 
tions in the two-level system are substantial and we can- 
not ascribe it a direct classical counterpart. As a conse- 
quence, we will demonstrate that the semi-classical pre- 
dictions are not recovered by introducing measurement 
induced decoherence of the oscillator subsystem (electro- 
magnetic field) in our model. We attempt to understand 
this lack of correspondence by showing that fundamen- 
tally non-classical characteristics, such as qubit-field en- 
tanglement, survive the decoherence. A weak coupling to 
an environment or a measurement device is conveniently 
described in terms of a master equation [19l [20] . In an 
unraveling approach to solving the master equation the 
Lindblad terms constitute imaginary state decay together 
with stochastic quantum jumps. Averaging over all these 
quantum trajectories results in the evolved state of the 
system. Single stochastic trajectories, which can be inter- 
preted as single experimental measurement realizations, 
are shown to support localization. As will be explained 
further below, it is important to distinguish this localiza- 
tion from dynamical localization [TJ JS] jSTJ [22] which is a 
quantum coherent effect that is destroyed by decoherence 
not generated by it £| . The localization phenomenon we 



discuss here only appears in the states of single stochastic 
trajectories and not in the full quantum probability dis- 
tribution (i.e. the average over an ensemble of measure- 
ments). We argue that signatures of chaotic behaviour 
in circuit QED can be seen in the photon statistics of the 
transmission line output field. 

This paper is organized as follows. The next section 
introduces the model system and the decoherence stem- 
ming from the position measurement. In section III we 



derive a set of semi-classical equations-of-motion where 
the field variable is treated at a mean-field level while 
a quantum treatment is maintained for the qubit. The 
results of the full quantum model are presented in sec- 
tion [IV| We especially focus on dynamics in phase space. 



including localization appearing due to ^selective mea- 
surements\ and methods to detect signatures of quantum 
chaos. Finally, we conclude in section [Vj 



II. PHYSICAL MODEL 

Our physical model consists of a superconducting qubit 
coupled to a driven transmission line resonator as, for 
example, in the experiments described in references [TT 
and [23 . The undriven setup is described by the Rabi 
Hamiltonian [24] (^ = 1) 



Hb = uj 






Q. 



CFz + V^gxcFj,, 



(1) 



where uj and ^ are the mode and qubit transition fre- 
quencies respectively, g the effective qubit-light inter- 
action strength, x and p the two field quadratures of 
the resonator mode, and cFx = |2)(1| + |1)(2| and cFz = 
1 2) (2 1 — |1)(1| are the Pauli matrices acting on the two 
qubit states |1) and |2). We will work with dimensionless 
variables where hw sets the energy scale, i.e. from now 
on uo = 1. The above Hamiltonian can be microscopi- 
cally derived from a minimal coupling model where the 
two- level, single mode, and dipole approximations are im- 
posed, and the "self-energy" of the field is neglected, see 
Ref. [25]. In the rotating wave approximation (RWA) one 
recovers the celebrated Jaynes-Cummings model 1 241126]. 
In traditional microwave and optical cavity QED setups, 
^ <C 1 such that this approximation is well justified [27] . 
while recent advances in circuit QED make it possible to 
reach regimes where the RWA breaks down [28 . When 
one goes beyond the RWA limit, especially for g > V^/2, 
the photon population no of the system ground state in- 
creases rapidly (no ^ g'^)- For the traditional configura- 
tion, namely a two-level system interacting with a single 
resonator mode, this build-up of photon population no in 
the ground state is mainly due to the fact that the self- 
energy of the Hamiltonian ([l]) has been left out [29] . This 
problem can be circumvented by following the timely 
suggestion made in reference [15]: by adding an exter- 
nal drive for the qubit one obtains an effective model 
of the Rabi form, but where the qubit-field coupling is 



scaled with the pump amphtude while the self-energy is 
unaffected, and, therefore, the two terms can be tuned 
independently. Hence, the self-energy can be neglected 
even in the ultrastrong coupling regime of the effective 
Rabi model. This fact is important for us here because 
we find that chaotic structures appear most definitively 
in the ultrastrong coupling regime. 

Making use of the Z2 parity symmetry of the Rabi 
Hamiltonian, the model was recently shown to be quasi- 
solvable over its entire parameter range [30 . This in- 
tegrability is somewhat surprising since the model does 
not possess a continuous symmetry. Moreover, integra- 
bility typically implies absence of chaos and the dynam- 
ics is therefore regular for all energies and parameters. 
However, as we will demonstrate, when the system is 
externally driven it may become chaotic for certain time- 
dependent drive amplitudes. The driven Rabi Hamilto- 
nian is given by 



HdR = a'^a + —az + ^ (a''' + a) a^ 



^r]{t) (ate-^^^^ + ae'^^^) , 



(2) 



where we have introduced the photon annihilation and 
creation operators a — {x^ip)/\/2 and a) = {x — ip)/\/2^ 
respectively. The drive amplitude will be taken to be 



r](t) = r]ocos{ujct), 



(3) 



i.e. the amplitude is modulated in time at angular fre- 
quency ujc between the maximum/minimum values of 
±7^0- The amplitude modulation frequency should not 
be confused with the central frequency Ud of the drive 
field upon which it imposes side bands. Returning to the 
quadrature representation, we have 



H, 



dR 



(^p^V2r]{t)sm{udt)y 



x^ 

y 



n 



{V2gc 



V2,ft)cos(c..t))f- ^^^^)-f^"^^\ 

(4) 
which will be our starting Hamiltonian. 

To gain further insight, it is convenient to perform an 
adiabatic diagonalization where the Hamiltonian is ex- 
pressed in the eigenvalue basis of the two-level matrix 
V{x) = Q&z/^ + gx&x with corresponding adiabatic po- 
tentials 



F±(x, ^) = y + V2r]{t) cos{uJdt)x ± ^ ^ 



n^ 



2g^x^. (5) 



The Born-Oppenheimer approximation consists of ne- 
glecting any coupling between the resulting adiabatic 
eigenstates [31 . If we let 7^0 = we note that the lower 
potential V^{x) has a double- well structure whenever 
g > VQ/2. The driving induces a back and forth rock- 
ing of the double-well. This type of driven double-well 
system mimics the Duffing oscillator [16 . However, for 



non-zero uud the driving also induces a "shaking" of the 
momentum. From Eq. Q it follows that this shaking 
effect can be accorded a time-dependent gauge potential 
A{t) = ^/2r]{t) sin{ujdt) giving rise to a synthetic electric 
field E = —dtA{t). While interesting, in this paper we 
shall not discuss this gauge aspect further. As will be 
described in the next section, the phenomena we are in- 
terested in stem from the rocking of the potential and this 
synthetic electric field will not change our arguments. 

Finally we want to include the impact of quadra- 
ture measurements [32]. In the Schrodinger picture 
and within the Markov approximation a standard ap- 
proach for describing decoherence is given by considering 
a master equation where the irreversible processes are 
attributed to a set of Lindblad operators j^j . Instead of 
considering field and qubit decay, we consider decoher- 
ence arising from the measurement of the x quadrature 
of the field. Thus, we work with a Markovian master 
equation of the form 



d 
dt 



p = i[p,HdR] ^ jC[x]p, 



(6) 



where p is the density operator for the full system, and 
the Lindblad operator 

C[x]p = K, (2xpx^ — x^ p — px'^) = —K, [f , [f , p] ] (7) 

accounts for the irreversible loss of coherence due to the 
measurement of the quadrature x [32] . The information 
loss rate is given by the parameter k,. Intuitively, a large 
K, is expected to produce a localization effect upon the 
quantum state. However, if the measurement back-action 
becomes too large it will inevitably affect the variable 
conjugate to the measured one. Thus, in order to achieve 
a localization effect, k should be in between the two ex- 
tremes. In the present bi-partite model this picture will, 
however, be very different. Note that the form of £[x] im- 
plies only loss of coherence and no dissipation. This will 
allow us to identify any localization effect in phase space 
as being due to decoherence rather than a reduction of 
the available phase space due to energy loss. 

A further remark before proceeding: whenever the sys- 
tem at hand is the subject of decoherence, information 
about it is irreversibly lost leading to an effective non- 
linear theory that can become chaotic in a classical sense, 
i.e. for short times nearby initial states move apart ex- 
ponentially with a rate set by the maximum Lyapunov 
exponent [34]. However, in an analysis like ours, any irre- 
versible process in quantum mechanics derives from the 
effective nature of the theory. Normally, it is assumed 
that the measurement device is macroscopic so that in- 
formation about the system is distributed over a large 
number of degrees of freedom and a reversed flow of in- 
formation back into the system is very unlikely. This 
justifies the Markovian assumption of infinitely fast de- 
cay of correlations within the measurement device [O] [20] . 
It is nevertheless an approximation and by waiting long 
enough information will flow back into the system. This 



waiting time will typically be very long in comparison to 
any experimentally relevant time scale so that the physics 
will appear to be non-linear in the experiment. 



III. SEMI-CLASSICAL ANALYSIS 

At large pumping, 770 > 1, we replace the mode vari- 
ables by their corresponding c-numbers as obtained in 
a coherent state ansatz. Then, the resulting effective 
Hamiltonian acts on the qubit sub-space alone, while its 
parameters depend on the boson field and must be de- 
termined self-consistently. In such a mean-field approach 
we neglect any quantum correlations between the qubit 
and the field, i.e. the qubit remains in a pure state and 
its presence does not imply decoherence on the large am- 
plitude resonator field. This allows us to derive a closed 
set of equations-of-motion, and thereby characterize the 
underlying semi-classical dynamics, i.e. whether it is reg- 
ular, chaotic, or mixed. 

Since the qubit is decoupled from the field, we can de- 
scribe it via two variables. After denoting the qubit state 
as Ix) = [/3ie-^^i , /32e-^^2]^^ ^^ere ^j and (j)j {j = 1, 2) 
are real, we introduce the (conjugate variables) inversion 
Z = P\ — ^2 and relative phase A^ = (/>i — (/>2- The semi- 
classical equations-of-motion become 



X = p + V2r]{t) sm{ujdt)^ 



p = —X 



gV2^yl- Z'^cos{A^) - V2r]{t) cos(cjrft). 



Z = gV2xVl^Z^M^^). 



n 
2 



- gV2x cos{A(i)) 



Vi-z^' 



(8) 
where dot denotes time-derivative. We have left out the 

decoherence induced Langevin forces [19 in this mean- 
field analysis. As a result, the semi-classical dynamics is 
independent of k which reflects the fact that the quadra- 
ture measurement only induces decoherence and no dis- 
sipation. From Eq. ([8|, we find the semi-classical Hamil- 
tonian 



Hsci{x,p,Z,A^) = 



p^ 



-^ pV2r]{t) sin{ujdt) 



-\-xV2r]{t) cos{uJdt) 



^, 



+gxV2Vl - Z2 cos(A0) + -Z. 

(9) 
We stress that the mean-field approximation is only ap- 
plied to the boson degree-of-freedom, in Eq. ([8| the spin 
is treated quantum mechanically. 

From Eq. ([8| we directly identify a trivial stable fixed 
point given by Xss = Pss = whenever r]o = (i.e. un- 
driven system). At g^ = r^/4, there is a Pitchfork bifur- 
cation where the fixed point turns unstable and two new 



stable fixed points arise {g'^ > r^/4), Xss = ^9\l^ — |^- 
These correspond to the minima of the double-well po- 
tential V~^{x). In an adiabatic picture, the driving in- 
duces an effective time-dependent "detuning" between 
the two double-well minima. For moderate driving, this 
dynamics can be treated as repeated Landau- Zener tran- 
sitions [35]. For large pumping rjQ (typically Ir^ol > 1^1), 
however, the double-well structure may be lost at the 
extreme of the driving, i.e. when \ri{t)\ attains its maxi- 
mum. In this case, the lower adiabatic potential oscillates 
between having one or two local minima. 
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Figure 1: Stroboscopic maps at t = T (a) and t = 4T (b). 
In (a), after one period T, the field phase space distribution 
still displays some structures, while in (b) they have more 
or less fully disappeared. The dimensionless parameters are 
Q = cjcz = 1, cjc = g = 1-5, and 770 = 3. The number of 
simulated trajectories Nd = 100 000. 

In this paper we use the semi-classical equations-of- 
motion to perform an analysis of the dynamics within the 
truncated Wigner approach (TWA) [37 . That is, we sim- 
ulate Nci ^ 1 trajectories of the above Eqs ([8| with ini- 
tial conditions for x and p randomly picked according to a 
Gaussian probability distribution corresponding to a co- 
herent state \a) with amplitude a = Xss/V^. The qubit 
initial state is taken A0(t = 0) = tt and Z{t = 0) = 0. 
This implies that the initial state predominantly popu- 
lates the right minima of the lower adiabatic potential 



V~^{x). The set of differential equations (p| are solved 
using the 4th order Runge-Kutta algorithm as modified 
by Gearys method^ which is especially suited to stiff equa- 
tions. 

Throughout this work we will consider the resonant 
situation Q = Ud = uj = 1, the pump amplitude r^o = 3, 
amplitude modulation frequency oOc = 1-5, and effective 
atom field coupling g = 1.5. This choice means that we 
are well into the ultrastrong coupling regime g > V^/2, 
and, as we see below, r^o = 3 is large enough to guarantee 
classical chaos. It should be pointed out that the results 
are not sensitive to the particular parameter values as 
long as the dynamics is chaotic. In fact, by fine tuning the 
drive parameters, chaotic evolution is also attainable in 
the "normal phase" g < ^/U/2 where the lower adiabatic 
potential has a global minimum at x = 0. 

With ujc = SuJdf^^ the Hamiltonian is periodic with 
period T = 47r. We consider stroboscopic maps of the 
phase space distribution of the resonator mode at times 
t = T and 4T. As can be seen in Fig. [l] after one pe- 
riod the distribution still shows some visible structures, 
while after four periods almost all such structures are 
washed out. For this choice of initial state and system 
parameters, the dynamics is completely chaotic lacking 
any islands of regular motion [34] . 

For a time-independent system, chaotic dynamics char- 
acterizes a smearing out of the distribution over the ac- 
cessible phase space energy shells. We note that the cor- 
responding effect in our time-dependent system occurs 
at rather short time-scales, e.g. the classical period of 
the harmonic oscillator Td = 2it in scaled dimensionless 
units. It is clear that the stroboscopic maps of Fig. [l] 
show different characteristics compared to those of the 
quantum Duffing oscillator [36]. This results from the 
additional degree-of-freedom of the spin. The effective 
phase space is here four dimensional instead of two di- 
mensional as it is for the Duffing oscillator. Therefore, 
the typical stretching and squeezing structure of the dis- 
tribution [36] is not as clear in the stroboscopic maps 
shown in Fig. [l] since they are projections onto the two 
dimensional xp-plane. Still, in Fig. fl] (a) some of the 
complex four dimensional phase space structure can be 
imagined. 



IV. QUANTUM ANALYSIS 

The master equation (|6| gives the evolution of the 
quantum state p under the inffuence of a certain form 
of decoherence which can be interpreted as a continuous 
position measurement. Extracting {x) implies measuring 
the field quadrature for an ensemble of equally prepared 
initial states [27 . The statistical average of such homo- 
dyne measurement sequences gives {x). Similarly, when 
solving (|6| numerically it is practical to transform the 
problem into one of solving a sequence of independent 
stochastic Schrodinger equations, each one of them re- 
sulting in a solution |'0i(t)) (note that this is a pure state). 



In this unraveling method any expectation value becomes 



►oo TV 



' Y.'^mt)m^{t)). Ev- 



{A)t = Tr [Ap{t)\ = lim^_ 

ery single {A)i = {tlji{t)\A\ipi{t)) can be viewed as one of 
the measurement sequences. As a consequence, the pho- 
ton detection is assumed to represent a perfect projective 
measurement where a 'click' in the detector is represented 
by a single 'quantum jump' in the stochastic Schrodinger 
equation. It is important to appreciate the distinction be- 
tween the outcome from a single measurement sequence 
and one deriving from an ensemble of measurements. We 
note that both are of experimental relevance. 

For the numerics we employed the split- operator 
method [38] combined with the quantum Monte Carlo 
method [39^. The split-operator method relies on fac- 
torizing, for short times St, the evolution operator into 
a real-space and a momentum part. The correspond- 
ing two propagators are applied in real and momentum 
space respectively. The method becomes exact in the 
limit 6t ^ 0. In general, the inherent time-scales set the 
constraints on the size of St. The final time is t = KSt, 
where K the number of steps used in the propagation. 
We varied 6t in order to check for convergence of our 
extracted results. The split-operator method allows us 
to simulate the coherent evolution, while for decoherence 
we apply the quantum Monte Carlo method, in which 
N different trajectories \i^i{t)) are simulated and where 
for each simulation i we periodically interrupt the coher- 
ent propagation and randomly choose either the absence 
or occurrence of a quantum jump characterizing a pho- 
ton detection [27 . This interruption is done K times, 
i.e. it agrees with the time-steps used for the coherent 
propagation, and we have verified that the jump prob- 
ability Pjump <^ 1 in each step. Note that a quantum 
jump in the present homodyne scheme is represented by 
multiplying the state \i^{t)) with x (and not with a as 
in a normal detection scheme where the photon inten- 
sity leaking out of the cavity is measured). As an initial 
state we will consider a coherent state predominantly oc- 
cupying one of the two minima of V~^{x,t = 0). Such a 
state can be approximated as |'^±(t = 0)) = | ± Q^)|0±), 
where a = Xqs/V^ is a real coherent state amplitude 
and |B±) = [±1 , l]^/>/2 is the initial qubit state. The 
±-sign represents the left or right potential minimum, 
respectively. This state approximates the ones used for 
the initialization of the TWA considered in the previous 
section. Notice that the initial state represents a disen- 
tangled field-qubit state. 

A Markovian master equation of the Lindblad form 
may imply (depending on the actual form of the Lind- 
blad terms) dissipation and decoherence. In a driven 
system, dissipation may balance the pumping such that 
a steady state containing a finite average number of pho- 
tons is approached [19 . The main result of decoherence, 
on the other hand, is to suppress off-diagonal terms of 
the density operator p [40 . This state "reduction" can 
hinder some evolution, similar to the quantum Zeno ef- 
fect. For chaotic systems, the role of reservoir induced 




quantum system Hi]. Moreover, the expectation values 
can become sensitive to small fluctuations in the initial 
states when decoherence is included, a typical charac- 
teristic for chaotic dynamics. At the same time, it has 
long been known that decoherence is the main source 
for loss of ' quant umness', e.g. quantum entanglement or 
negativity of the Wigner function [9^, "20^ . It might there- 
fore not come as a surprise that decoherence might help 
in obtaining a proper quantum-classical correspondence 
for chaotic systems. While the issue of decoherence in 
chaotic systems has been rather well studied for cases in 
which the system has an obvious classical limit, for sys- 
tem with no clear classical limit little is known ^18^ . The 
situation is even more complex in bi- or multi-partite sys- 
tems where entanglement may persist between the differ- 
ent subsystems even in the presence of decoherence [42j . 



Figure 2: (Color online) The Wigner phase space distribu- 



tion (10) for the resonator field and for four different decay 



rates n ^ 0, 0.005, 0.01, 0.05 (a)-(d) respectively The effect 
of decoherence is seen to have several consequences; i over- 
all smoothening - washing out sub-Planck structures, ii slight 
suppressing of phase space spreading, and Hi abate negativ- 
ity. The rest of the dimensionless parameters are as in Fig. [l] 
with a propagation time t = 200 for all four plots. 




Figure 3: (Color online) Identical to Fig.[2]but displaying the 
Husimi distribution (11) instead. In contrast to the Wigner 



distribution, the Husimi distribution is positive definite and 
lacks sub-Planck scale structures due to the Gaussian averag- 
ing as discussed in the main text. 



decoherence has been discussed in numerous works, see 
for example [2l [36l [41]. Decoherence stemming from 
quantum measurements has been explored as well, see 
Ref. [171 132] • In the case of the quantum Duffing oscilla- 
tor, it has been verified that the expectation values {A)t 
obtained from the corresponding classical model agreed 
much better with those of the quantum model when 
decoherence is included than those given by the closed 



1. Phase space analysis 

Quantum phase space distributions are a practical tool 
for achieving a deeper intuition for quantum dynamical 
systems [19, 43 . The idea is to visualize the quantum 
evolution in a phase space picture. The Wigner distribu- 
tion for a quantum state p{t) (in our case the resonator 
boson mode) is defined as [43] 



W{x,p,t) = ^J 



dy{x-y/2\p{t)\x + y/2)e'Py. (10) 



Because it does not have to be positive everywhere in 
phase space, Wix^p^t) is not a proper probability dis- 
tribution. Typically, negativity of W{x^p^t) is identified 
with non-classical phenomena such as quantum coherent 
superpositions. The Wigner distribution is normed to 
unity, JJdxdpW{x,p,t) = N_{t) + A^+(t) = 1, where 
N±{t) measures the negative (— ) and positive fraction 
(+) of the distribution. Note that normalization does 
not imply \N±{t)\ < 1. The marginal distributions of 
W(x,p, t) give the real and momentum space probabil- 
ity distributions, i.e. they are positive definite. Another 
characteristic of W(x,p, t) is that it allows for sub-Planck 
structures [44 . As a counterpart to classical phase space 
distributions, the Husimi Q-function seems more as a 
natural choice. The Q-function is positive definite and 
in contrast to the Glauber P- functions, it does not be- 
come singular. It is defined as [191 US] 



Q{x,p,t) = -{a\p{t)\a), 



(11) 



where |a) is a coherent state with amplitude a = ar -\- 
iai = V2{x + ip) {ar and a^ are hence both real). It 
follows that the Husimi Q-function is in fact a Gaussian 
averaged Wigner function [45], 

Q{x,p,t) = - f f dxdpe-^''-^^'-^P-P^'w{x,p,t). (12) 



This smoothening not only extinguishes the negativ- 
ity, but also prevents sub-Planck structures from form- 
ing [46^. Interestingly, despite this averaging, the Q- 
function still contains all the information there is to know 
about the state p{t). While the Q-function has the ap- 
pealing property of being positive definite, its marginal 
distributions do not reproduce the correct real and mo- 
mentum space probability distributions. Note that the 
phase space distributions are defined for general states 
p(t), with no restriction on purity. This applies espe- 
cially to us where we are interested in the distributions 
for the field. The density operator for the field pfit) 
is typically in a mixed state even without coupling to 
the environment since it is obtained from the full density 
operator by a trace over the qubit's degrees-of- freedom. 

For a chaotic system like the present one, we expect 
the phase space distributions to show an ergodic collapse 
in which they spread out over the accessible phase space. 
Even at resonant driving, cj^^ = (^ = 1, the closed model 
has a bounded energy shell due to the oscillating pump 
amplitude rjit) and the coupling to the qubit, which is 
something also indicated by Fig. [l] After the collapse, a 
seemingly irregular phase space structure builds up. In 
a closed quantum system, the above behaviour is typical 
for quantum thermalization [47 . How the phase space 
structures develop in a chaotic open system is, however, 
not well explored. In our model we cannot talk about 
quantum thermalization since the system is both driven 
and exposed to decoherence. 

We have numerically calculated W(x,p, t) and 
Q{x,p,t) for different decay rates hz. The results are 
shown in Figs. [2] and [3] respectively. Comparing the two 
figures, it is apparent that the sub-Plank structures in 
the Wigner distribution {h = 1) do not occur in the Q- 
function. Furthermore, as expected decoherence enters 
as smoothening of the Wigner distribution, and for large 
Hi the sub-Planck structures are completely lost in the 
Wigner distribution (Fig. [2](d)). 

Instead of deriving the distributions directly from the 
state pf{t)^ it can alternatively be obtained from the 
corresponding Fokker-Planck equations [I9j. This ap- 
proach gives us additional insight into the phase space 
dynamics. The full bi-partite Q-function [19], Q{z^a) = 
TT^fi+Ul^) (^7 ^\pi'^)\^i ^)i where \z) is a spin coherent state 
with complex amplitude z ^48j, obeys the equation 

dtQ{z, Of, t) = (£ + /:qdiff + >Cmdiff) Q{Z, Of, t) . (13) 

The first differential operator on the r.h.s represents a 
classical drift 



ir]{t)e-'^^'da + ida 



a-\-3g 



z -\- z 



1 



-^id, [-nz + g{l - z^){a + a*)] + c.c, 
the second quantum diffusion 

>Cqdiff = igdadz{l - Z^) + C.C, 



(14) 



(15) 




Figure 4: (Color online) Evolution of the effective phase space 
widths /\xp (a) and 5xp (b) for various decay rates. The black 
line in the upper plot is obtained from the semi-classical TWA 
result. The remaining lines (solid red, dot-dashed blue, dotted 
green, and dashed purple) display the full quantum results. 
The width /\xp of the full phase space distribution shows a 
weak Ai:-dependence. In contrast, single trajectory widths 5xp 
has a strong Ai:-dependence which demonstrates a localization 
effect. The rest of the dimensionless parameters are the same 
as in Fig. [l] 



and the final term describes the measurement induced 
diffusion 



£r 



idiff 



l^{da-da'f ■ 



(16) 



From Q{z^a^t) we find the field Q-function by integrating 
out the spin coherent state variable z. The form of the 
second term £qdiff prevents sub-Planck structures from 
forming in agreement with the Heisenberg uncertainty 
relation [49 , i.e. the origin of this term is the quantum 
pressure which becomes very large whenever the distri- 
bution is tightly squeezed. This explains the quantum 
nature of this effect (and this term indeed vanishes in 
the classical limit h ^ 0). Note further that >Cqdiff mixes 
spin and field variables. The last term £mdiff stems from 
the noise resulting from the position measurement and 
only contains field variables. In this model, where x is 
being measured, it has a diffusive effect on its conjugate 
variable p. Its form is the same in the Fokker-Planck 
equation obeyed by the Wigner distribution which ex- 
plains the smoothening effect seen in Fig. [2] as a^: is in- 
creased. To summarize: C reproduces the classical dy- 
namics, >Cqdiff gives quantum corrections, e.g. prevents 
sub-Planck structures from forming and generates qubit- 
field entanglement, and finally i2mdiff describes the back- 
action of the measurement. 



As already mentioned, chaotic evolution typically im- 
plies a rapid spreading of the phase space distribution. 
The 'size' of the distribution can be estimated from the 
width 



A^^ = AxAp = V((x2) 



{m{{p')-{p)'). (17) 



i.e. the product of the quadrature uncertainties. Figs. |2] 
and |3] give the distributions as evolved according 
to the master equation ([6|, but as discussed above 
we may also consider the distributions Qi{x^p,t) and 
Wi{x,p,t) for single trajectories i. Each one of 
these, corresponding to single stochastic simulations of 
the unraveled master equation, gives a width S[ 



(i) 



the average width 



(p)f) which are used to define 



1 



N 



i=l 



(18) 



Axp and dxp represent, respectively, the non-selective 
width in which measurement results are not recorded and 
selective width in which the measurements are projective 
and information about the results of a measurement are 
fed back into the system. For the same decay rates tz as 
in Figs. I2] and Isl we display the time-evolution of A^p 
and Sxp in Fig.^ (a) and (b). For a comparison to the 
semi-classical results we have in addition calculated Axp 
within the TWA framework using ([s]) . The figure reveals 
numerous interesting phenomena: 

1. Even for large decoherence rates hz, the position 
measurement does not localize the state Pf{t) such 
that it can be said to approximate the semi-classical 
result. This is in stark contrast with results found 
for the Duffing oscillator [32]. In Ref. [18] a simi- 
lar breakdown of the quantum-classical correspon- 
dence was conjectured to be due to the quant um- 
ness of one of the bi-partite subsystems (two-level 
system). Below we will motivate this by calcu- 
lating the qubit-field entanglement, and it will be 
shown that despite strong decoherence substantial 
non-classical correlations between the qubit and the 
field survives. 

2. A much more distinct localization effect is found in 
the selective width 6xp- For large times, the quan- 
tum distributions can even remain more localized 
than the semi-classical one given by the TWA. We 
have verified that the localization occurs in both 
conjugate variables x and p, and not only in the 
measured position x. 

The localization here occurs due to quantum 
jumps. Every time a jump occurs, the wave func- 
tion tp{x^t) is multiplied by x. Such a 'projective 
measurement' interrupts the evolution and these 
stochastic 'kicks' of the wave function cause a de- 
structive interference leading to the localization ef- 
fect. This is reminiscent of dynamical localization. 



which is the analogue in time of ID Anderson local- 
ization in space, and occurs because of random (but 
coherent) interference that prevents the wave func- 
tion from spreading (the latter has in fact been seen 
in a periodically "kicked" optical cavity filled with a 
non-linear medium [51] ) . Here the kicks result from 
projective measurements. However, it is interesting 
to note that the wave function of a single simulated 
trajectory ilji{x^t) is a pure state and so the mech- 
anism for dynamical localization is perhaps quite 
close even if the origins are different. Meanwhile, 
the full quantum probability distribution, as given 
by the solution of the unraveled master equation, 
and which also coincides with the average over all 
trajectories, does not show the same localization 
effect. This is not surprising since different trajec- 
tories are added incoherently and so one should not 
expect destructive interference. 

3. In both Fig. [2] (a) and (b), we find that the ef- 
fect of the decoherence appears after some time td - 
This dwell time is a characteristic of open chaotic 
systems [53 . The dwell time characterizes a pe- 
riod during which the open system is transparent 
to decoherence and evolves similarly as the corre- 
sponding closed system. From Fig. [4] (a) we see 
that td ^ 20 irrespective of hi. Normally, during 
the dwell time the system evolves essentially classi- 
cally. In our case we again find a breakdown from 
the typical behaviour as the semi-classical broad- 
ening diverges from the quantum one long before 
the dwell time is over. 

One explanation for the failure of the semi-classical de- 
scription to reproduce the fully quantum ones even in the 
presence of decoherence comes from the fact that some 
quantum properties are maintained despite the decoher- 
ence. Since the state p{t) is in general mixed, standard 
entanglement measures like the von Neumann entropy 
cannot be applied to estimate qubit-field entanglement. 
Neither can the concurrence be utilized since the two sub- 
systems do not constitute two qubits. As a measure of 
entanglement we use instead the negativity [54 



En{p) = -Y1 



\\i\-\i 



(19) 



where A^ is the ith eigenvalue of the partial transpose 
p^^ of p [55]. The definition (19) means that the nega- 
tivity is the sum of negative eigenvalues A^. In principle, 
En{p) < is only a necessary (but not sufficient) condi- 
tion for entanglement, but, nevertheless, we will take it 
as an indicator of quantum correlations. We have a lower 
bound En{p) > 1/2. The negativity for a relatively large 
K, and for early times is displayed in Fig. p](a). At later 
times, En{p) shows rapid variations around the value 
^ —0.05. Thus, the entanglement is not completely de- 
stroyed by the decoherence. It has been shown that for 
closed quantum systems, non-zero entanglement can be 
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Figure 5: (Color online) The evolution of the negativity ^Ar(t) 
(a) and negative fraction N-(t) (b) for k = (black solid 
lines) and k, = 0.05 (red dashed lines). For larger times, the 
negativity follows a similar behavior with declining oscillation 
amplitudes, i.e. the bi-partite system remains entangled even 
for long times despite the non-zero k. Decoherence suppresses 
entanglement, but a substantial fraction survives even large 
decay losses. Similar behaviour is encountered for the neg- 
ative fraction of the Wigner distribution as demonstrated in 
the lower plot. However, it seems that N-{t) is more sensitive 
to decoherence. The remaining parameters are as in Fig. [l] 



a signature of quantum chaos [6l[56]. In general, chaotic 
dynamics shows large entanglement in closed systems. 

Typically, given a Lindblad term jC[A] with [A, ^] = 
the state p{t) tends to relax to a steady state which 
is the ground state of A^A [57 . Here, since [x^H] ^ 
there is an interplay between the two terms, unitary 
Hamiltonian time-evolution and non-unitary decoherent 
time-evolution, and the state is not expected to approach 
some pure steady state. Indeed, the state p{t) is highly 
mixed which can be verified by calculating the purity 
P{t) = Tr[p^(t)] [58 . Although not shown here, we have 
calculated P{t) for different hz's and only in the special 
case where z^: = is the state pure, i.e. P{t) = 1 for /^ = 
and P{t) < 1 otherwise. 

As mentioned above, another signature of non-classical 
features is the negative ratio N- {t) of the Wigner func- 
tion |59! . The rule-of-thumb is that decoherence tends to 
suppress non-classical effects, such as; squeezing, sub- 
Poissonian photon statistics, and coherent superposi- 
tions. We also saw an example of this in Fig. [2] Exam- 
ples of the negative fraction A^_ (t) is presented in Fig. [s] 
(b). While N-{t) is suppressed, there is always a fraction 
that survives. Compared to the negativity, N_{t) seems 
somewhat more sensitive to decoherence. 



2. Detection 

We have seen examples of how the effects of deco- 
herence come into play for chaotic quantum dynamics. 
Therefore, a most relevant question concerns experimen- 
tal observation of both quantum chaotic dynamics it- 
self, and the effect decoherence has upon it. This far 
we have considered the phase space distributions and es- 
pecially their widths. The full phase space distribution 
is in principle attainable experimentally via state tomog- 
raphy [60, £lj, see also [6 . Such homodyne detection 
relies on experimental repeatability, i.e. repeated mea- 
surements on copies of equally prepared systems. Since 
the system is manifestly chaotic and non-linear, small 
fluctuations in the preparation will cause large varia- 
tions in the evolved state and the measured results may 
seem random. The measurement result will therefore 
more likely represent a state averaged over many "re- 
alizations" . This, however, is in itself a signature of the 
chaotic nature. Thus, state tomography should render 
a phase space distribution with seemingly random struc- 
tures. 

We may ask whether there are other quantities that 
are more easily obtained and still show chaotic charac- 
teristics? It turns out that even the statistical properties 
of the intensity of the output cavity field contain infor- 
mation of the quantum chaotic behaviour. The phase 
space distribution, apart from spreading, will bunch back 
and forth along the x-direction. When passing the cen- 
ter at x = 0, population transfer between the two qubit 
components takes place, i.e. adiabaticity breaks down. 
For short times, when the wave packet is still localized, 
this will be manifested as sudden jumps in (a^). Af- 
ter the delocalization, the clear transitions between the 
two components turn into rapid fluctuations. This will 
also be reflected in large fluctuations in the field intensity 
(a^a). Classically, the ergodic evolution implies that a 
phase space trajectory frequently intersects all accessible 
"rings" with radii r = y^p^ + x^/2 in phase space. Thus, 
a classical solution spends some time at all accessible dis- 
tances from the origin of phase space, and such chaotic 
motion results in pronounced field intensity fluctuations. 
For a coherently driven system, the intensity of the out- 
put field at time t, (aout^out)^, leaking the cavity mirror 
is proportional to the cavity mode intensity {a^a)t [62] . 
Consequently, direct photon detection of the output field 
gives a handle on the cavity photon fluctuations. 

In Fig. [6] we compare the photon variance of the full 
field state pf{t) (i.e. not for single trajectories) for a con- 
stantly pumped system (a), and with the one driven with 
an oscillating amplitude (b). In (a) we show the variance 



^n = V ^(^^^) ^ ~ (a^a)'^ and in (b) the scaled variance 
An/{d^a). The reason why we do not divide the variance 
in the regular case (a) is because (a^a) periodically be- 
comes approximately zero and would cause divergences. 
For the constant driving case we find a nicely oscillating 
behaviour, while in the chaotic case the variance shows 



10 



much greater fluctuations. Another difference is that de- 
coherence enhances fleld intensity fluctuations in the reg- 
ular regime, while it suppresses it in the chaotic regime. 
As a comparison we also display the TWA result. As is 
seen, without decoherence the quantum variance exceeds 
the semi-classical one. For strong enough decoherence, 
the quantum variance may actually be smaller than the 
semi-classical variance. 
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Figure 6: (Color online) Variance of the photon intensity. The 
upper plot (a) gives the variance for a system pumped with a 
constant amplitude, i.e. the dynamics is regular. The lower 
figure (b) shows instead the scaled variance for the system 
with an oscillating drive. In both plots the red dashed line 
gives the full quantum result with no decoherence, and the 
green solid line gives the same except that k = 0.01. The black 
solid line displays the corresponding TWA result. Comparing 
(a) and (b) clearly indicates a more irregular structure in 
the photon variance whenever the system is evolving in the 
chaotic regime. The rest of the parameters are as in Fig. [l] 



demonstrated that the set of semi-classical equations-of- 
motion predicts classical chaos. In general, decoherence, 
as it is known to diminish quantum coherence, may ex- 
plain how classical chaos can appear from quantum me- 
chanics. This is typically the case for systems with a 
classical limit. However, we showed that this method of 
achieving classical-quantum correspondence fails in our 
bi-partite model where one of the two subsystems is man- 
ifestly quantum. In particular, within a wide range of de- 
coherence rates k we found that the semi-classical results 
could not be retrieved. By calculating the negativity we 
concluded that this discrepancy between quantum and 
semi-classical chaotic evolution could originate from the 
presence of truly quantum properties such as entangle- 
ment even in the presence of strong decoherence. While 
the full quantum state p{t) did not show much localiza- 
tion, the states appearing from 'selective measurements' 
gave a form of localization distinct in origin from stan- 
dard dynamical localization but actually rather close in 
mechanism to the stochastic Schrodinger equation we em- 
ployed to describe the measurements. We further argued 
that the chaotic nature of the system should show up in 
the photon statistics which are readily measured in exper- 
iments. We note that the amplitude of the information 
transfer rate hz should be controllable by adjusting the 
homodyne detection. 

It should be pointed out that our results apply to other 
chaotic models exposed to decoherence. Maybe the most 
relevant example is the the Dicke model [63 describing 
a set of N two-level qubits coherently interacting with 
a single cavity mode. The ultrastrong coupling regime 
was recently achieved experimentally in an effective Dicke 
model [64 , and by exciting the system sufficiently it can 
be expected to enter a chaotic regime, where our predic- 
tions should be detectable in the output cavity field. Our 
analysis assumes pure decoherence and no dissipation, 
but we have also tried other forms of the Lindblad terms 
including, for example, photon decay, and found similar 
results. When the system is dissipative, one must keep 
in mind that photon decay itself may cause a localization 
of the phase space distribution due to energetic consid- 
erations that are different from the localization discussed 
here. 



V. CONCLUSION 



In this paper we have proposed using a circuit QED 
setup for studies of quantum chaotic dynamics in a sys- 
tem with no clear classical limit. Special attention was 
paid to the infiuence of decoherence deriving from con- 
tinuous quadrature measurements of the resonator field. 
Since the measurement back-action only induces decoher- 
ence and no dissipation, the corresponding semi-classical 
equations-of- motion are transparent to this effect. We 
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